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ABSTRACT 



We report observations of SN 2004dt obtained with the Very Large Telescope 
of the European Southern Observatory on August 13.30, 2004 when the supernova 
was more than a week before optical maximum. SN 2004dt showed strong lines of 
O I, Mg II, Si II, and Ca II with typical velocities of absorption minimum around 
17,000 km s _1 . The line profiles show material moving at velocities as high as 
25,000 km s _1 in these lines. The observations also reveal absorption lines from 
S II and Si III with a velocity of only 11,000 km s _1 . The highest velocity in the 
S II features can be traced to about 15,000 km s _1 , much lower than those of O, 
Mg, Si, and Ca. SN 2004dt has a polarization spectrum unlike any previously 
observed. The variation of the polarization across some Si II lines approaches 2%, 
making SN 2004dt the most highly polarized SN la ever observed. In contrast, 
the strong line of O I at 777.4 nm shows little or no polarization signature. The 
degree of polarization points to a silicon layer with substantial departure from 
spherical symmetry. A geometry that would account for the observations is one 
in which the distribution of oxygen is essentially spherically symmetric, but with 
protrusions of intermediate-mass elements within the oxygen rich region. 
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Subject headings: supernovae - cosmology - spectropolarimetry 



1. Introduction 

As for core-collapse supernovae, high-quality spectropolarimetry of Type la supernovae 
(SNe la) is revealing a complex array of behavior that has opened a new window of explo- 
ration of this class of explosions (Wang et al. 1996; Wang, Wheeler, Hoflich 1997; Howell 
et al. 2001; Wang et al. 2003a; Wang et al. 2004). While there are indications that the 
degree of polarization of core-collapse supernovae increases with time after explosion (Wang 
et al. 2001, 2002; Leonard et al. 2001), the time evolution of Type la supernovae seems to 
be rather different (Wang et al. 2003a). 

Polarimetry of SN la shows that the degree of polarization of normal SN la is < 0.3%, 
which is much lower than the typical polarization of core-collapse supernovae (Wang et al. 
1996). Spectropolarimetry of the normally-bright SN 1996X at optical maximum tentatively 
showed intrinsic polarization of about 0.3% which implies geometrical asphericities of about 
10% (Wang et al. 1997). Hydrodynamic instabilities just beyond the layers of complete 
thermonuclear equilibrium may have contributed to the polarization signal from SN 1996X. 
Howell et al (2001) showed that SN 1999by, a sub-luminous SN la, was polarized at a level 
of 0.7% in the red continuum. 

The observations of the normally-bright SN 200 lei revealed the geometric structure of 
an SN la in different chemical layers. The degree of polarization showed significant time evo- 
lution (Wang et al. 2003a). An exceedingly strong Ca II feature at velocity of 20,000 - 26,000 
km s _1 was highly polarized. This feature showed that the outermost layers of the supernova 
ejecta had a highly aspherical geometry (Wang et al. 2003a; Kasen et al. 2003). Deeper in- 
side the ejecta, the silicon-rich layer also displayed asphericity as evidenced by the detection 
of polarization of ~ 0.3% for the Si II 635.5 nm line. When the photosphere receded deep 
inside the iron-rich layers, the polarization diminished and the supernova appeared spher- 
ical. Wang et al. (2003a) listed several possible origins of this asymmetric high-velocity 
material: matter from a disk or companion star that is swept up by the supernova ejecta; 
newly-synthesized material produced by thermonuclear reactions in the high velocity collision 
between supernova ejecta and the surroundings; and a clump of matter ejected from the in- 
trinsic thermonuclear combustion process. Gerardy et al. (2004) examined the high-velocity 



*Based on observations collected at the European Southern Observatory, Chile (ESO Progr. No. 073. D- 
0565(A)). 
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Ca II feature in SN 2003du (for which there was no spectropolarimetry) and concluded that 
the strength of the Ca II in that event was consistent with solar abundance of calcium in 
a hydrogen plasma (that the bulk of the matter was helium or carbon/oxygen cannot be 
ruled out) that was swept up from distances consistent with the expected separation of the 
companion star. 

We present in this paper observations of SN 2004dt. This event showed exceptionally 
strong and broad lines of Si II, and O I and revealed the strongest polarization signal we have 
yet measured for a SN la. The discovery and pre-maximum observations are described in §2, 
the spectral and spectropolarimetry features and their implications for ejecta geometry are 
discussed in §3. The inferred ejecta geometry is given in §4. The observations of SN 2004dt 
are compared to other SN la in §5. Comments on theoretical models are given in §6, based 
on the observational results of SN 2004dt. We give in §7 some discussions on the impact of 
asymmetry on the use of SN la as cosmological probes. Finally, in §8, we provide a brief 
summary and some general discussions of the most important conclusions we have derived 
from spectropolarimetry. 



2. Discovery and observations 
2.1. Discovery 

SN 2004dt in NGC 799 was discovered on an unfiltered image taken on August 11.48 at 
a magnitude of about 16.1 (Moore & Li 2004). Salvo, Schmidt & Wood (2004) reported that 
a spectrum taken on August 12.72 UT showed SN 2004dt to be a young SN la a few days 
before maximum light. The photospheric expansion velocity deduced from the absorption 
minimum of Si II 635.5 nm was about 16,500 km s _1 , adopting the recession velocity of 5915 
km s^ 1 reported by the NED for the host galaxy NGC 799. Patat, Pignata, and Benetti 
(2004) reported that a spectrum obtained on Aug. 13.17 UT showed some unusual features. 
They noted that the O I 777.4 nm transition was particularly pronounced, with a broad 
symmetric P-Cygni profile and also had an expansion velocity of about 16,800 km s -1 . They 
also report that the Ca II H & K line gave 22,200 km s _1 (we think this feature is mostly 
due to a Si II line, see §3.2). Patat, Pignata, and Benetti (2004) suggested that the high 
velocities of the lines indicated a very early phase. They noted that other intense lines visible 
in the spectrum were from Fe III and Si III. 
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2.2. Observations and data reduction 

We observed SN 2004dt with the ESO-VLT on 2004 August 13.30 UT with the FORS1 
instrument in polarization mode. From preliminary light curves of Li et al. (2004, private 
communication), this corresponds to 6-8 days before optical maximum. Four exposures of 
1200 second durations were taken with the waveplate at position angles of 0, 45, 22.5, and 
67.5 degrees. Flux and polarization standard stars were observed to calibrate the flux levels 
and polarization position angles. The data were reduced in a way similar to that described 
in Wang et al. (2003a). 

2.2.1. Spectrophotometry 

The observations consist of multiple exposures with the airmass varying from 1.02 to 
about 1.20. The slit size was 1 arcsec and was aligned with the paralactic angle. The seeing 
was below 1 arcsec which allows for approximate spectrophotometry calibrations. Fig. 1 (top 
panel) shows the total flux spectrum. Starting in the blue, there is no sharp drop of flux blue 
ward of the strong Ca II H & K absorption as commonly observed in normal SN la at optical 
maximum. The Ca II H & K feature appears very strong although, as we argue in §3.2, it 
should probably be identified as Si II 385.9 nm. The strong line of Si II 635.5 nm shows 
the typical blue-shifted absorption that characterizes SN la at 600 nm, but it is especially 
strong, broad, and fast moving. The minimum corresponds to a velocity of about 17,000 
km s -1 , and the blue wing extends to about 25,000 km s" 1 . We also observe a strong line of 
O I 777.4 nm as pointed out by Patat, Pignata, and Benetti (2004). This line also showed an 
absorption minimum at about 17,000 km s^ 1 with the blue wing again extending to about 
25,000 km s -1 . To first order, the Si and O occupy the same velocity space. Our data 
show two components of the Ca II IR triplet, a broad absorption at 820 nm and a "notch" 
indicating a separate, high-velocity component at about 800 nm. These represent velocities 
of ~ 16,000 and 23,000 km s _1 , respectively. The lower velocity component corresponds 
approximately to the velocity minimum of Si II and O I. The higher velocity component 
seems kinematically distinct as in SN 2001el (Wang et al. 2003a) and SN 2003du (Gerardy 
et al. 2004) and to correspond roughly to the blue wing of the Si II and O I lines. Despite 
the likely two-component nature, the highest velocity, again, can be traced to 25,000 km s _1 
for the Ca II IR triplet, similar to the other high- velocity lines. 

The strong absorption feature at 421 nm is due to Mg II 448.1 nm, blended with Fe III 
and Si III lines. If this feature is predominantly Mg II, then the absorption minimum of the 
Mg II line is blueshifted by 17,000 km s _1 , consistent with the lines of Si II, Ca II, and O II. 
The strength of the 421 nm feature is comparable to the O I 777.4 nm line, as measured by 
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the ratio of line minimum to local continuum. 

The absorption minima of the S II "W" lines are measured to be at 523.5 and 543.1 nm 
which implies a velocity of only 11,000 km s" 1 . This is significantly lower than the velocity 
of the Si II lines. Weak lines at 439.5, 454.2, and 464.0 nm, tentatively identified with Si III 
456.0 nm, S II 471.6 nm, and S II 481.6 nm, respectively, also show the much lower velocity 
of 11,000 km s -1 . The presence of Si III 456.0 nm at 11,000 km s _1 , not at 17,000 km s" 1 
shows that doubly ionized Si exists only close to the central region. Silicon is the only element 
present in two different ionization stages. They share the emerging bimodal distribution of 
radial velocities: Singly ionized silicon has a bulk velocity of 17,000 km s _1 whereas Si III 
expands with 11,000 km s _1 . 

Table 1 gives a schematic overview of the kinematic properties of the main ions observed. 

2.2.2. Spectropolarimetry 

The second, third, and fourth panels of Fig. 1 show the polarization, P, and the reduced 
Stokes vectors, U, and Q. The degree of polarization is corrected for biases due to observa- 
tional noise using the equations derived in Wang, Wheeler, & Hoflich (1997). The Q vector, 
in particular, shows a total excursion of ~ 2 percent with sharp wavelength dependence, 
indicating strong, intrinsic polarization. The bottom panel of Fig. 1 gives the polarization 
position angle. It is nearly constant from 330 nm to 630 nm, but then shows a strong change 
with possibly significant excursions at wavelengths longer than 630 nm. This change could 
be due to a combination of polarization by interstellar dust and the intrinsic polarization by 
the supernova ejecta. Note that the polarization position angle is periodic with respect to 
the position angle on the sky with a period of 180°. 

Fig. 2 gives the data in the Q/U plane where each point represents the Stokes parameters 
rebinned to 25 A intervals. The bin size was chosen to be larger than the spectral resolution 
to minimize correlations among neighboring lines. The total excursion is nearly 3 percent. 
The scatter of points is predominantly due to line features. The data points span a range 
of about 1.5 percent in Q and more than 2 percent in U. Fig. 2 shows that the bulk of the 
points fall along the line that represents the "dominant axis" of the geometry (Wang et al. 
2001), determined by making a principle component decomposition of the data points on 
the Q/U plane. The scatter orthogonal to this axis suggests that there are smaller but real 
deviations in the geometry from this dominant axis that may also constrain the physics of 
the explosion. 
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3. The distinguishing observational features 

3.1. The spectral lines and profiles 

The colors synthesized from the spectrum give U — B ~ -0.6. This blue color (as verified 
by the strong blue flux in Fig. 1), suggests that this event, if approximately "normal," is 
close to or before U maximum. This may mean that at this phase, SN 2004dt has not yet 
been subject to the strong Fe II absorption leading to the characteristic UV deficit of normal 
SN la near maximum light (Wheeler & Harkness 1992; Kirshner et al. 1993). 

The strong lines of Si II 635.5 nm and O I 777.4 nm in the total flux spectrum suggest 
that there is substantial ejected silicon and oxygen, moving at velocities approaching 25,000 
km s _1 . This is consistent with a significant amount of incompletely burned matter. The lack 
of a UV deficit suggests that there is not much iron in the photosphere at the phase of these 
observations. To investigate further the possibility of an iron deficit, we have performed line 
identifications using the SYNOW code (Fisher et al. 1999). Assuming a composition with 
only Si II, S II, O I, and Mg II, we found that the observed spectral features are reasonably 
well explained without requiring iron. What distinguishes SN 2004dt from other Type la 
supernovae appears to be that a larger than usual amount of silicon is required to explain 
the deep spectral features from Si II. The SYNOW line identification is shown in Fig. 3. 
The absence of Fe features is another striking property of the spectrum of SN 2004dt at this 
epoch. 

The presence of lines at a velocity of 11,000 km s _1 suggests that even at this early 
phase, the ejecta is not Thompson opaque at this lower velocity. The velocity of the Si III 
455.3, 456.8, and 457.5 blend is only 11,000 km s -1 , implying that doubly-ionized silicon 
exists mainly in this lower- velocity region. The low velocity S II line suggests that S is more 
concentrated in the low- velocity region than Ca, O, Mg, and Si. 

3.2. The interstellar polarization toward SN 2004dt 

From the strong modulation of the raw polarization shown in Fig. 1 it is immediately 
evident that the intrinsic linear polarization of SN 2004dt is exceptionally large. The vectorial 
nature of polarization can lead to confusing conclusions if the interstellar polarization (ISP) is 
not properly determined and subtracted; unfortunately, the determination of ISP inevitably 
requires some model assumptions. We have tried to avoid assumptions that are tied to the 
models of SN la explosions. 

The central assumption made here is that the ISP dominates where the continuum 
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polarization of the supernova is the lowest. The effective continuum polarization probably 
gets decreased by overlapping spectral lines. Since there are many more such lines in the 
blue portion of the spectrum than in the red portion (Wang et al. 2001; Howell et al. 2001; 
Wang et al. 2003a), we have selected the regions marked in Fig. 1 to represent the minimal 
intrinsic polarization; they avoid significant discrete absorption features. 

It is seen in Fig. 2 that the deduced ISP corresponds to 0.2 ±0.1% in Q and —0.2 ±0.1% 
in U. This low value is consistent with the absence of an interstellar Na I D line (§3.6) if the 
gas-to-dust ratio of the host galaxy is normal. 

In Fig. 2, the data points used for the determination of the ISP are embedded in a 
large cloud of other low-polarization wavelength bins. If, following Wang et al. (2001), a 
principle-component decomposition technique is used to define the dominant and the or- 
thogonal polarization axes, the ISP will, therefore, be close to the dominant axes. As Wang 
et al. (2001) have shown, some minor error in the determination of the ISP does not grossly 
distort a qualitative interpretation if the true ISP is close to the dominant axis. 

A similar constellation of data points prevailed in SN 2001el (Wang et al. 2003a). The 
key corroborating observation was that at later phases the polarization of the continuum as 
well as of spectral features retreated along the dominant axis and, starting a week after max- 
imum, remained within a similar cloud of low-polarization wavelength bins. A preliminary 
analysis of later-epoch data of SN 2004dt shows a similar behavior. 

The resulting ISP-corrected decomposition into dominant and orthogonal axis is pre- 
sented in Fig. 4. 

3.3. The intrinsic polarization of SN 2004dt 

Ultimately, the choice of interstellar polarization must depend on the totality of obser- 
vations and a self-consistent picture of the physics of the explosion and radiative transfer 
through the complex ejecta. Important details of the interpretation of the data may yet 
depend on the particular placement of the ISP. For this paper we have attempted to restrict 
our discussion to features that are relatively independent of the choice of the ISP. 

SN 2004dt shows a polarization spectrum unlike any SN la yet observed. Compared 
to the well-observed SN 200 lei, the Si II lines are especially prominent in the polarization 
spectrum, but the lines of the Ca IR triplet are less so. 

With the above choice of the interstellar polarization, we see from Fig. 4 that the 
absorption minima of the P-Cygni profiles correspond in general to the peaks of the spec- 
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tropolarization profiles. Along the dominant axis polarization spectrum, the polarization in 
the range 620 - 720 nm is around -0.3% for the choice of interstellar polarization adopted in 
Fig. 4. This wavelength range shows no strong spectral lines and it is likely that continuum 
polarization by Thompson scattering dominates this wavelength region. 

The high quality data we obtained allow for a detailed comparison of the spectral and 
spectropolarimetry profiles of several distinct lines, as shown in Fig. 5. Spectropolarimetry 
allows us to revisit the issue of line identification as discussed already in §3.1. 

3.3.1. Si II lines 

One of the strongest polarization features in Fig. 4 and Fig. 5 is the spike at 600 
nm which corresponds to the Si II 635.5 nm line. The strength of this feature, ~ 2%, is 
unprecendented in the observed polarization spectra of SN la. There is also a strong spike 
in the polarization at 365.0 nm, the location of the "Ca H & K" line. We suspect this is 
due primarily to Si II 385.0 nm that is only slightly contaminated by the Ca lines. Note 
that the amplitude of this polarization feature, ~ 2%, is nearly identical to that of the Si II 
635.5 nm line. No other features in the polarization spectrum show this level of polarization. 
The spectral and spectropolarimetry profiles of Si IIAA385.9, 635.5 nm lines are shown in 
Fig. 5(a) and (c), respectively. The polarizations are the highest at the absorption minima 
of these two features. The spectral and spectropolarization profiles of the two lines are 
remarkably similar. The peaks of polarization occur at a blueshift of about 20,000 km s _1 for 
both lines then drop sharply at higher velocity. This strengthens our identification that they 
are from the same ionic species. The dominant-axis spectrum shows another feature at 390 
nm that we attribute to Si II 413.0 nm. Note, however, that the same feature was tentatively 
identified as C II 426.7 nm with an expansion velocity of about 24,100 km s _1 by Patat et 
al. (2004). Our identification with Si II would put the line at the same velocity of about 
17,000 km s" 1 as the other Si II lines. 

The line at 480 nm can be due to the blend of Si II 504.1, 505.6 nm, and Fe II 491.3, 
501.8, 516.9 nm. As can be seen in Fig. 4, the contributions from Fe II must be small as 
evidenced by the absence of the strongest Fe II line in this wavelength region, Fe II 516.9 
nm, in the total flux spectrum. The small feature at 570 nm is likely to be due to Si II 595.8, 
597.9 nm. To the red of the Si II 635.5 nm line, the polarization dips and is then constant 
to about 720 nm. With these identifications, all the Si II lines have velocities from 15,000 
to 17,000 km s _1 . 

The lines of Si II 413 nm, Mg II 447.1 nm, and Si II 504.1, 505.6 nm all have polarization 
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~ 1% in Fig. 4, significantly less than the strongest polarized lines of Si II, but suggesting 
some commonality in their geometry. All these lines have similar dominant polarization 
position angles. 

3.3.2. Si III and S II lines 

Weaker lines from Si III, and S II are also polarized, at a level of about 0.3%, and 
with polarization position angles similar to those of the Si II lines. The blueshifts of the 
absorption minima for these lines are around 11,000 km s _1 which is significantly lower than 
the strong Si II lines. 

3.3.3. Mg II 420.0 nm 

The polarization feature at 420 nm is mostly due to Mg II 447.1 nm. The profiles are 
shown in Fig. 5(b) for this line. Like the polarizations of the Si II 635.5 nm and Si II 385.9 
nm lines, the polarization of the Mg II 447.1 nm line peaks also at a blueshift of around 
20,000 km s _1 , but the polarization does not exhibit a sharp drop at higher velocities and 
appears to extend to even higher velocities than those of the Si II lines. 

3.3.4. Call lines 

The polarized spectral feature at 820 nm is clearly due to the Ca II IR triplet; the 
absorption minimum of the Ca II feature shows a velocity of 16,000 km s -1 , consistent with 
that of the Si II lines. Note that, unlike SN 2001el which showed a highly polarizing detached 
Ca II shell at a velocity of 22,000 km s _1 , the high- velocity component of the Ca II IR triplet 
of SN 2004dt with a similar velocity is not particularly strong nor does it have an especially 
high polarization. SN 200 lei and SN 2004dt are distinctly different in this regard. The lower 
velocity component of the Ca II IR triplet in SN 2004dt shows a polarization spike of about 
0.7-2%, comparable to the Mg II line. The higher velocity component shows a spike of about 
0.3%, comparable to the S II and Si III lines. These polarized components of Ca II seem to 
share approximately the same polarization angle as all the other polarized lines, ~ 150°. 

At the same velocity as the Ca II IR triplet, only a weak dip is observed for the 
Ca II H&K lines. The weak Ca II H&K lines do not produce a prominent polarized feature. 
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3.3.5. I 777.4 nm line 

Despite the very strong spectral line of O I 777.4 nm in the total flux spectrum, there is 
little sign of that feature in the polarization spectrum. The spectral and spectropolarimetry 
profiles of this line are shown in Fig. 5(c). This low polarization is in sharp contrasts to all 
the other strong lines from Mg II, Si II, and Ca II. 

3.4. Summary of observed features 

The rather diverse spectral information extracted from our data can be summarized 
as follows (see also Table 1): (1) Lines of Mg II, O I, Ca II, and Si II are at the highest 
velocities, their absorption minima are typically blue shifted by 17,000 km s _1 , and the 
lines can be traced to a velocity of around 25,000-28,000 km s -1 ; (2) The S II lines have 
absorption minima at only 11,000 km s" 1 ; the highest velocity of these features is less than 
15,000 km s -1 ; (3) The absorption minimum of the Si III 455.3, 456.8, and 457.5 blend is 
blue shifted by 11,000 km s" 1 ; (4) All the high- velocity features except that from O I 777.4 
nm are significantly polarized, at levels from 0.7% to 2.2% with two features of Si II being 
especially distinct; (5) Weak lines from S II, and Si III are also polarized, at a level of about 
0.3%, but share the same polarization position angles as those of the other polarized lines; 
(6) the high- velocity shell of Ca II has a velocity comparable to the blue wings of the strong 
lines, but with a polarization amplitude and angle comparable to S II, and Si III. 

4. The Inferred Ejecta Geometry 

4.1. Bubbles and Plumes 

A critical issue is whether the observed polarization indicates a clumpy ejecta or a 
globally aspherical structure with little in the way of small scale features. The Q-U diagram 
provides a powerful tool to study the dumpiness of the ejecta. If the asymmetry is caused by 
an axially symmetric structure, the resulting polarization will enscribe a straight line on the 
Q-U diagram. This is easy to understand as the polarization position angle is only related to 
the position angle of the symmetry axis on the sky. Deviations from the straight line on the 
Q-U diagram are expected only if the structure departs from axial symmetry, for instance 
by being fragmented and clumpy. 



The observed polarization for each line can be fit by a linear relation: 

U = a + (3Q. 



(1) 
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Such a fit gives the dominant axis of that particular line in the Q-U plane. The goodness 
of this fit is a direct measure of the smoothness of the ejecta. For the strongest lines, the 
region for this linear fit was chosen to be from -25,000 km s _1 to -10,000 km s~ x for Si II 
635. 5nm, 385. 9nm, and Mg II448.1 nm; this corresponds closely to the absorption minima 
of the P-Cygni profiles. The corresponding range was chosen to be from -25,000 km s _1 to 
-13,000 km s -1 for the O I 777.4 nm line to avoid contamination by telluric features. The 
situation is less clear for weaker lines which are often blended with other lines. For the S II 
"W" feature, the velocity range was taken to be from -15,500 km s _1 to -9,500 km s _1 and 
both lines were used for the linear fit. For Ca II IR triplet, the lines are blended and the 
range of the fit was taken to be 786.0 nm to 834.2 nm. 

Fig. 6 shows the Q-U diagrams for the strongest lines, Si II 385.9 nm, 635.5 nm, 
Mg II448.1 nm, and O I 777.4 nm. Table 2 gives the details of the linear fits for all the 
lines that allow for a linear fit on the Q-U plane. The columns in Table 2 are: the atomic 
line, the blueshift range of the data used in the fit, a, (3, the resulting position angle of the 
straight line, the x 2 , ,he degrees of freedom of the fit, and the probability that the association 
with a straight line is by chance. 

The linear fit provides a good description to the observed Q-U diagram only for the 
Mg II and the S II lines. For the other ionic lines, the fits are highly inconsistent with a 
straight line. It is interesting to note that the O I line, despite its lack of a well-defined 
axis, is inconsistent with a random Gaussian distribution and therefore the data actually 
suggest intrinsic polarization of the oxygen, but with no well-defined symmetry axis. The 
Si II and Ca II lines all show clear evidence of a dominant axis, but with significant intrinsic 
deviations from the dominant axis. Only the Mg II gives a good fit (x 2 per degree of freedom 
~ 1) indicating a distribution along a single, well-defined axis with no measureable deviation 
due to fragmentation or clumping. The Si II 635.5 nm line gives a fit with y 2 per degree 
of freedom 4, and a position angle that only differs from that of Mg II by about 2a. This 
suggests that the silicon may share the same dominant axis as the magnesium, but with 
significant clumping. 

This analysis suggests that the distribution of silicon follows some large scale asymmetry 
that gives rise to the dominant axis, but that there are also prominent smaller scale structures 
in silicon that lead to significant deviations from the dominant axis in the Q-U diagrams. 
The polarization of the Si III line at lower velocity but with the same position angle suggests 
that this asymmetric silicon structure extends from at least the photosphere at about 11,000 
km s" 1 all the way to the highest observable velocity for the Si II lines, i. e., 25,000 
km s -1 . The transverse scale of this silicon-rich structure should be comparable to the 
size of the ejecta where the Thompson optical depth is around 1, so that the asymmetric 
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line opacity of the structure will have a large enough area covering factor to produce both 
the strong absorption and the large degree of polarization. This radially-extended silicon 
structure is fragmented or clumpy. The size of the small scale structures that give rise to 
the departure from the dominant axis must be larger than 1,200 km s _1 , which corresponds 
to the approximate resolution of the data, and to the bin-size of the data points on the Q-U 
diagrams given in Fig. 6. Our data are not sensitive to clumps smaller than 1,200 km s _1 . 

The fact that the Mg II line is consistent with a straight-line fit in the Q-U plane with 
no significant deviation from the principal axis means that, although its geometry follows 
quite nearly that of the Si II globally, the magnesium has a smooth large structure with 
no noticeable small scale clumps or fragments. This has several implications. First, the 
success of the straight line fit for this line implies that the errors associated with the Q and 
U vectors are reasonable; this lends confidence in the robustness of the fits for the other 
lines. Secondly, this implies that the distribution of magnesium is less clumpy than that of 
silicon. This could be understood if some amount of magnesium coexists with the silicon-rich 
clumps but that there are magnesium-rich regions that do not contain noticeable amount of 
silicon. The number density contrast of magnesium in different clumpy regions (presumably 
generated by explosive carbon burning) must be small; the number density contrast of silicon 
in those different clumpy regions must be large. 

The lack of a clear dominant axis for the O I line and the obvious inconsistency of the 
data with random Gaussian distribution on the Q-U plane suggests that the distribution 
of oxygen is even more pervasive than magnesium, but the opacity must be very clumpy, 
perhaps due to large variations of oxygen abundances in different chemical regions. 

The Q-U diagram shows undoubtably that at the highest velocity, the ejecta of SN 2004dt 
are highly asymmetric, clumpy and fragmented. 

4.2. Nucleosynthesis at the Highest Velocities 

The null detection of polarization across the O I 777.4 nm line demands further expla- 
nation. From the spectral profile, we know that oxygen co-exists with silicon, magnesium 
and calcium. Why does not the O I 777.4 nm line share the same polarization properties of 
lines from the other elements? 

A simple explanation is that oxygen is distributed in a different way from silicon, mag- 
nesium and calcium, although they share the same velocity space. Oxygen is likely more 
pervasive in the outer region, but the protrusions of products from nuclear burning are dis- 
tributed in holes or filaments in this region. Such protrusions also contain oxygen. In such a 
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geometry, the polarization from oxygen and from the other elements are drastically different 
as their geometrical structures are not identical but occupy different regions. 

If so, then the null detection of polarization across the O I 777.4 nm line implies that the 
oxygen rich region does not exhibit the same asymmetric geometry as the magnesium/silicon 
rich regions. This would require that the magnesium/silicon rich regions also contain a 
significant amount of oxygen so that oxygen appears to be more "spherically" distributed. 
Indeed, in explosive carbon burning, oxygen is expected to co-exist with magnesium, a 
principal burning product. 

The observed high degree of polarization of high- velocity magnesium suggests that the 
magnesium-rich structure was created by explosive carbon burning. The fact that the mag- 
nesium structure is polarized with about the same amplitude and position angle as the silicon 
lines also shows that magnesium is not confined to magnesium-rich, silicon-poor regions, but 
also exists in silicon-rich regions. This in turn suggests that the oxygen that must be present 
outside the magnesium/silicon structure in order to have a very small net polarization must 
be the pre-explosion oxygen from the progenitor white dwarf. 

The unusual polarization of SN 2004dt is closely related to the observed velocity struc- 
tures of the various ions. In the outermost portions, pre-explosion oxygen co-exists with 
newly-synthesized Si, Mg, and Ca. This outermost layer moves at a velocity around 17,000 
km s" 1 and the highest velocity approaches 25,000 km s _1 . The silicon-rich layer extends 
deep inwards to a velocity of around 11,000 km s -1 , as evidenced by the presence of Si III 
and S II lines at that velocity. The broadness and unusually high blueshift of the Si II lines 
are evidence that these lines are optically thick. Iron and sulfur are noticeably absent in the 
high-velocity layer. The strong polarization observed across the P-Cygni lines of the newly- 
synthesized elements is the result of the asymmetric distribution of the chemical elements 
in this this layer. The polarization position angles of these lines are similar, suggesting that 
they share a common geometrical structure; but with oxygen permeating the whole struc- 
ture, magnesium occupying a globally asymmetric region, and silicon occupying a similar 
globally asymmetric region but, in addition, being broken into clumps or fragments. 

The structure of the ejecta is sufficiently asymmetric that spectral lines can be formed 
in significantly different velocity regions. The attribution of the polarization to asymmetric 
chemical distributions requires the structures to have considerable radial extent because the 
angle of the dominant axis is the same for lines formed at different depths, that is, for both 
Si II at 17,000 km s _1 and Si III at 11,000 km s" 1 . The photosphere in the conventional sense 
of Thompson optical thickness probably has a complicated geometry that will complicate 
attempts to understand the radiative processes. For simplicity, one may still assume a 
photosphere that is located at a velocity much smaller than the 17,000 km s _1 as derived 
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from the absorption minima of the P-cygni lines of Si II lines. The weaker S II "W" lines 
and Si III lines at 439.5 454.0, and 464.0 nm may provide a hint of the location of the 
photosphere. If so, at least along certain lines of sight the photosphere is at a velocity as low 
as 11,000 km s _1 . This rather low photospheric velocity may not be as surprising. From light 
curve fitting, we know that the pre-maximum luminosity of SN la rises approximately as t 2 
(Goldhaber et al. 2001), which implies a photosphere at a fixed velocity if the atmosphere 
is approximated by a blackbody. A photospheric velocity of 11,000 km s _1 agrees with the 
observed luminosity and light curve behavior of typical SN la. 



4.3. Continuum polarization and the geometry of the photosphere 

We note that most of the discussion in this paper relies on polarization variations across 
spectral features, not on the level of continuum polarization. We will not elaborate on the 
geometrical structure of the photosphere. With the complicated geometry shown by the 
various lines, the continuum may be formed in a complicated clumpy environment. It is 
likely not to have a smooth geometry like the one assumed for ellipsoidal models. The 
fact that the emission peaks of the strongest P-Cygni profiles are all blueshifted by about 
4,000 km s -1 (see Fig. reflineprofies) suggests that there is a substantial amount of electron 
scattering at high velocities. The rather extended region of electron scattering effectively 
blocks the receding side of the ejecta and can shift the emission peak of the P-Cygni profile 
to the blue. This interpretation of the strong blueshift has to be reconciled with the presence 
of faint spectral line features at velocities as low as 11,000 km s" 1 through more detailed line 
profile modeling. Lacking such models, it should be remarked, however, that a clumpy ejecta 
is able to produce P-Cygni profiles with a blueshifted emission peak together with absorption 
features at low velocity (Wang & Hu 1994). In such a scenerio, the patchy opacity structure 
allows formation of spectral features at low velocity together with shielding of the receeding 
side of the supernova atmosphere. 

Taking the adopted components of interstellar polarization at face value, we can see 
that the level of continuum polarization can be as low as 0.2-0.3%, consistent with what was 
found for SN 2001el. The amplitude of continuum polarization is larger in the red part of the 
spectrum than in the blue part, but with the sign of the polarization flipped. The relatively 
low level of continuum polarization does suggest that the asphericity of the photosphere is 
not very large, perhaps only at a level of less than 10 percent as has been derived for other 
Type la supernovae such as SN 1996X (Wang, Hoflich, Wheeler 1997) and SN 2001el (Wang 
et al. 2003a). The exact level of continuum polarization awaits a full analysis of the data 
set which covers post-maximum phases of the supernovae. 
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5. SN 2004dt as a member of a sub-class of SN la 

5.1. Spectroscopically similar events 

SN 1997bp and SN 2002bo displayed high velocity absorption features similar to those 
of SN 2004dt. Strong polarized features were also observed across the Si II 635.5 nm line 
for both SN 1997bp and SN 2002bo. Discussion of the polarization features of SN 1997bp 
and SN 2002bo will be given in a forthcoming paper. We discuss here only their distinct 
spectroscopic behavior. As pointed in Benetti et al. (2004), supernovae such as SN 1984A 
(e.g. Branch 1987) and SN 2002bo show unusually high velocity lines. A comparison of the 
spectra of SN 1997bp, SN 2002bo, SN 2004dt, and some spectroscopically normal SN la is 
shown in Fig. 7. The characteristic behavior of SN 1997bp, SN 2002bo, and SN 2004dt is that 
they appear spectroscopically normal, but with the lines much broader and at much higher 
velocities. Benetti et al. (2004) suggest that for this group of supernovae the burning to 
silicon penetrated to higher layers than in more normal SN la. The high velocity, broad lines 
can be understood by the presence of a fast moving shell which is rich in intermediate mass 
elements and pre-explosion oxygen. The extension of this shell in velocity space requires 
detailed spectral modeling. 

Both SN 1997bp, and SN 2002bo show strong narrow Na I D absorption that implies 
that they suffer considerable amount of dust extinction from the host galaxy. The equivalent 
widths of the Na I D absorption are 1.5 A, and 2.5 A for SN 1997bp and SN 2002bo, 
respectively. No strong Na I D line is detected for SN 2004dt. For SN 1984A, Barbon et al. 
(1989) derived E(B-V) = 0.45±0.05. Kimeridze & Tsvetkov (1986) found A B = 1.2±0.2. 
Although these values are uncertain, the observed color around maximum does suggest large 
interstellar extinction. For SN 1997bp, Riess et al. (1998) give Ay = 0.62, suggesting 
E(B — V) ~ 0.2. A significant amount of dust along the light of sight of SN 1997bp in 
the host galaxy is also consistent with the presence of a strong Na I D line. Benetti et al. 
(2004) found E(B-V) = 0.43 ±0.10 for SN 2002bo, in agreement with the observed strong 
Na I D in SN 2002bo. SN 2004dt does not show a noticeable Na I D line, and final light 
curves are not yet available to estimate the dust content on the line of sight; however, our 
flux-calibrated spectrum shows B — V ~ 0.07 at around day -7. Although uncertain, this 
implies an E{B — V) of 0.14-0.2 when compared to colors of typical SN la at comparable 
epoch. Typical B — V for an unextinguished supernova at this phase is around -0.07 {e.g. 
Phillips et al. 1999). These values of E(B — V) are significantly larger than typical values 
of E(B — V) found for Type la supernovae (Knop et al. 2003). 

It is not clear to us whether the extinction of these supernovae is caused by circumstellar 
or interstellar dust. High resolution spectroscopy may be helpful in this regard, but it may 
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be difficult to use high resolution spectroscopy alone to tell whether the dust is within a 
few parsecs or a few hundred parsecs from the progenitor. Late time light echo observations 
(Patat 2004) combined with late time polarimetry may set more direct constraints on the 
location of the absorbing dust. A larger sample of SN 2004dt-like objects will resolve the 
issue of how strongly they are associated with dusty environments. 

5.2. SN 20024dt and other polarized supernovae 

SN 2004dt continues to add to the data base that strongly suggests that SN la are 
polarized in a very interesting fashion, but with a variation from event to event that we are 
just beginning to probe. 

The observations have revealed both high and low velocity spectral features about a 
week before optical maximum. Strong variation of polarization across spectral lines of inter- 
mediate mass elements, especially silicon, is observed; but most puzzling is non-detection of 
polarization in the O I 777.4 nm line. 

SN 2004dt is different than SN 2001el. For the latter event, the strongest polarization 
feature was the high-velocity component of the Ca II IR triplet with Si II polarized at a 
much lower level than the Ca II IR triplet. For SN 2004dt, the polarization of the Ca II 
IR features is prominent, but the polarization levels of the lines of Si II are the strongest. 
SN 2001el showed continuum polarization at the level of about 0.3%. SN 1996X tentatively 
showed similar levels of continuum polarization. As we discussed before, the precise level of 
continuum polarization of SN 2004dt is difficult to determine because of the highly-polarized 
spectral features. With the adopted level of interstellar polarization, the continuum level 
may be comparable to that of SN 2001el and SN 1996X; however, the large polarization 
across most of the spectral lines makes SN 2004dt distinctively different from SN 200 lei. 
In Wang et al. (2003a), we argued that it is important to study the correlations among 
the degree of polarization of the high- velocity Ca II IR triplet, the strength of the high- 
velocity component of the Ca II IR triplet, and the velocity of the Ca II IR triplet. The 
observations of SN 2004dt are consistent with a positive correlation between the strength 
of the high-velocity Ca II IR triplet and the degree of polarization. The observations also 
seem to confirm the speculation that higher velocity events should show larger polarization. 
Obviously, more spectral polarization data are needed. 
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5.3. SN 2004dt and Sub-Types of SN la 

All of the four Type la supernovae identified to have abnormally high velocity Si II 
lines show also evidence of unusually large reddening by the dust in the host galaxy. Both 
SN 1997bp and SN 2002bo are found to be significantly polarized, at levels around 1-2%. 
There were no polarimetry data on SN 1984A. This group of objects is likely to be char- 
acterized by high velocity, high dust extinction, and intrinsically large polarization before 
optical maximum. 

We do not yet know whether these high velocity Type la supernovae form a separate 
sub-group of SN la, or whether they are simply a familiar sub-group of objects viewed at a 
special angle. For the latter possibility, SN 1991T-like events may be candidates to be the 
counterparts viewed at different angles. The only other sub-group of SN la that is known 
to be associated with dusty environments are SN 1991T-like events. Perhaps SN 1991T-like 
objects are viewed in a direction with very little intermediate mass elements in the outer 
layers, whereas SN 2004dt-like objects are the same objects viewed in the direction with 
substantial pre-explosion elements on the outside. 

6. Comments on theoretical models of Type la supernova explosions 

SN la involve the combustion of a degenerate C/O white dwarf of the Chandrasekhar 
mass in a binary system (Whelan & Iben 1973; Nomoto 1980, 1982). Major theoretical 
questions concern the ignition and propagation of the flame through the white dwarf. Three 
major scenarios have been proposed: detonation (Arnett 1969; Hansen & Wheeler 1969), 
deflagration (Nomoto, Sugimoto, & Neo 1976; Nomoto, Thielemann, & Yokoi 1984), and 
delayed detonations (Khokhlov 1991; Yamaoka et al. 1992; Woosley & Weaver 1994). Three- 
dimensional calculations show that pure deflagration models leave unburned carbon and 
oxygen mixed into the center (Gamezo, Khokhlov & Oran 2002; Reinecke, Hillebrandt, 
Niemeyer, 2002). Delayed-detonation models leave very little unburned carbon and tend 
to produce the observed stratified composition structure (Gamezo, Khokhlov & Oran 2004; 
Marion et al. 2003; Hdflich et al. 2001, 1995, Wheeler et al. 1998). 

6.1. Deflagration 

In SN 2004dt, we see a highly asymmetric and clumpy structure in silicon that extends 
from at least as deep as 11,000 km s _1 to at least 17,000 km s -1 . This fragmented structure 
could, in principle, be consistent with the turbulent structure produced in 3-D simulations 
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of deflagration models (Gamezo et al. 2002; Ropke, Niemeyer, J. C, & Hillebrandt 2003). 
The layer at velocity around 17,000 km s -1 shows pre-explosion oxygen and fragments of 
intermediate mass elements, but with no iron. Deeper inside, sulfur is found at a velocity 
around 11,000 km s _1 which shares the same asymmetry signature as the high velocity, 
oxygen-rich layer. 

In deflagration models, silicon co-exists with iron; the absence of strong iron features 
in the spectrum may be inconsistent with published pure deflagration models, but we cau- 
tion that the exact amount of iron that can be hidden in the spectra without producing 
pronounced spectral features is uncertain and quantitative spectral modeling is needed to 
set the constraints on deflagration models. The velocities of the silicon and magnesium are 
significantly higher than predicted by pure deflagration models. It thus seems that turbulent 
deflagration models are capable of producing the observed clumpy composition structure, but 
are perhaps incapable of producing the kinematics of the ejecta. Deflagration models also 
predict unburned matter, C and O, mixed down to low velocities. Subsequent observations 
of this event will help to put constraints on this possibility. Deflagration models also tend 
to leave a layer of completely unburned C and O on the outside. It is very important to put 
limits on the amount of carbon in the spectra. 



6.2. Delayed-detonation 

An alternative is the delayed-detonation (DD) model (Khokhlov 1991). This model 
assumes that at a certain stage of deflagration, the nuclear burning becomes supersonic and 
therefore switches to a detonation. In order to have a high velocity partially burned layer, 
there must be significant pre-expansion (due to the subsonic deflagration stage) before the 
onset of detonation, which sufficiently reduces the density of the progenitor white dwarf 
before the nuclear burning flame propagates to the surface layers. It may also be possible 
that there is a layer on the surface of the progenitor white dwarf with a density that is too 
low to burn at all. Because of the low density near the surface, the products of the nuclear 
reactions are mostly intermediate mass elements. A strong detonation wave is unlikely to 
generate a turbulent silicon layer such as we apparently observe in SN 2004dt. One possible 
way out is that the detonation wave weakens significantly near the surface of the star, in 
which case the cellular structure associated with detonation instability (Boisseau, et al. 1996; 
Gamezo et al. 1999) can grow to significant scales. 

Alternatively, the detonation could be switched on at a stage where the deflagration 
flame has propagated sufficiently close to the low density part of the expanding white dwarf 
that the products of detonation bear the imprint of the complex structure of deflagration. 
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A delayed detonation would be able to boost the outer layers to the velocities we observe. 
It may also be consistent with the absence of iron features at very high velocity. Detailed 
modeling is needed in order to understand the clumpy structure of the products of nuclear 
burning, especially silicon, that this data indicates. 

6.3. Off-centered delayed-detonation 

A different model which naturally produces a globally aspherical Type la is given in 
Livne (1999) and more recently in 3-D calculations by Gamezo et al. (2004) in which the 
detonation is triggered at a point off the center of the white dwarf rather than in a spherical 
shell. The result is that the inner nickel region is offset from the centroid of the density 
distribution. The resulting off-center production of gamma-ray energy input can induce 
polarization even in a spherical density distribution (Chugai 1992; Hoflich 1995). In addition, 
the silicon is predicted to be left with a distribution that is displaced from the centroid of the 
density distribution. This leads naturally to the possibility of a globally asymmetric silicon 
distribution that has the potential to produce significant polarization by large scale blocking 
of the photosphere as well as smaller scale clumps of silicon resulting from the turbulent 
burning. This sort of model thus has the capability of producing all the polarized signatures 
of the silicon-rich layer. Since the photosphere might also be significantly aspherical from 
the energy input or the asymmetric density distribution, a large continuum polarization may 
be produced by such a model; this should be confronted with detailed theoretical modeling. 
This model can be tested by spectroscopic observations in the nebular phase. We will defer 
a detailed discussion of this model to a later paper when the analysis of the complete data 
set of our observations is undertaken. 



6.4. Gravitationally-confined detonation 

In the model recently proposed by Plewa, Calder & Lamb (2004), the explosion begins as 
a deflagration born slightly off-center in the white dwarf. Buoyancy dominates the evolution 
of the deflagration, resulting in a very rapidly rising plume of burned material (see also 
Livne, Asida & Hoflich 2005). This plume displaces the outer (unburned) layers of the 
star, which propagate ahead of burned bubble material across the surface of the star. This 
flow converges at the point opposite from where the plume emerges, and the resulting rapid 
compression of the unburned material triggers a detonation. An advantage of this model is 
that it may not require triggering the detonation artificially. This model has about 1 percent 
of a solar mass of iron group elements on the surface of the white dwarf resulting from the 
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surface flow of the plume material. As for pure deflagration models, detailed spectroscopic 
modeling is needed to tell if this amount of iron is allowed by observations. In addition, the 
requirement of significant pre-expansion so that sufficient intermediate mass elements are 
produced in the velocity range from 11,000 km s _1 to 25,000 km s _1 to match the observed 
spectra may further constrain this model. Note that pre-expansion is necessary to produce 
intermediate mass elements in the outer layers, but can significantly weaken the shock of the 
merging plume matter; this will work against the triggering of the detonation. Therefore 
these observations will set limits on the parameter space where this mechanism may work. 
Due to pre-expansion at the stage of detonation, it is also possible that this mechanism 
will produce highly aspherical chemical structures such that one side (where detonation is 
triggered) consists of mainly iron rich material on the outside, whereas the other side (where 
the burned bubble emerges) consists of mainly of intermediate mass elements. Depending 
on which side faces the observer, the target might appear as an iron-rich SN 1991T-like, 
or a silicon-rich SN 2004dt-like; however, the details of these constraints await thorough 
theoretical models of these data. 



7. Implications for Supernova Cosmology 

This study illustrates that the nature of Type la supernovae is more heterogeneous 
than is easily revealed by observing methods that are less telescope time intensive than 
spectropolarimetry. We have shown that polarimetry provides information that is often 
strongly complimentary to other observables and so can provide important discrimints among 
theoretical models of SN la explosions. 

In order to make SNe la better understood standard candles it will be essential (i) 
to obtain enough data to extract the core commonalities from all the individual cases and 
(ii) to find ways of inferring certain special conditions also from low-resolution spectra and 
photometry because only they permit SNe to be observed at cosmological distances. If 
unrecognized sub-types of SNe occur in different proportions at different redshifts or in 
different environments, supernova cosmology experiments could be severely compromised by 
such uncertainties. 

There are two effects of asymmetry on the observed magnitudes of a supernova. First, 
asymmetry affects the total integrated light from a supernova in a particular filter, therefore 
causing a view-angle dependence of distances. Secondly, asymmetry will cause intrinsic color 
dispersions which can make the determination of color-excess due to dust in the host galaxy 
difficult. The second effect can become very important as errors on color are amplified when 
extinction corrections are to be applied. 
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The nature of events like SN 2004dt may require this group to be separated from other 
supernovae, and a re-calibration of the methods (e.g. Phillips et al. 1999; Wang et al. 2003b) 
commonly used to derive distances and dust reddening of SN la. This supernova may repre- 
sent a particular sub-set of SN la, along with SN 1984A, SN 1997bp, and SN 2002bo, with 
especially high-velocity matter for which the explosion mechanism is significantly different 
than that for their lower velocity counterparts such as SN 1994D. 

An association with dusty environments would suggest also that these events are likely 
to be from younger progenitor stars, and that they may be more abundant at higher red 
shifts than in our local universe. Identifying these events and applying the proper relations 
for distance calibration is also critical for removing the effects of progenitor evolution. 

8. Summary and Discussion 

Spectropolarimetry is a straight forward, powerful diagnostic tool for supernova explo- 
sions. Contrary to conventional perceptions that the interpretation of polarimetry data is 
complicated and heavily model dependent, many important insights on the geometric struc- 
ture of supernovae ejecta can be derived from diagnostic analyses that do not rely on detailed 
modeling. 

In this study, we have presented high quality spectropolarimetry of the Type la SN 2004dt. 
These data allow us to probe the geometric structure of the ejecta from around 11,000 km s _1 
all the way to about 25,000 km s -1 . 

SN 2004dt adds to the examples of how these simple diagnostic tools can be applied to 
well observed supernovae to study the geometric structure of the ejecta. Spectroscopically, 
we found that most of the spectral features can be identified with intermediate mass elements 
with little or no contributions from iron. From spectropolarimetry, we found the high velocity 
regions of SN 2004dt to be highly aspherical, clumpy and probably turbulent. 

The results of this study provide significant constraints on theoretical models of Type 
la explosions. Polarimetry of SN la probes the origins of the intrinsic diversity of SN la, and 
is important for improving them as precise tools of cosmology. 
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Table 1: Summary of Observed Polarization Features 



Ion 


Rest wavelength(s) [nm] 


Radial velocity [km s x ] 


Polarization 






bulk 


blue edge 


degree [%] 


angle 


Cont 


350-600 


11,000 


N/A 


0.0 


N/A 


Cont 


600-850 


11,000 


N/A 


-0.3 


60° 


SII 


471.6, 481.6, 543.3, 545.9 


11,000 


15,000 


0.3 


150° 


Si III 


455.3/456.8/457.5 


11,000 




0.3 


150° 


Ca II 


850.0/854.2/866.2 


16,000 


27,500 


0.7-2 


150° 


Mg II 


447.1 


17,000 


28,000 


1 


150° 


01 


777.4 


17,000 


27,500 


<0.2 


N/A 


Sill 


385.0, 413.0, 504.1, 505.6 635.5 


15,000-17,000 


25,000 


0.5-2 


150° 
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Table 2: Properties of Polarized Lines 



Line 


V range (km S _1 ) 


a 


P 


P.A. b 


x 2 


DoF 


Prob. 


Si II 635.5 


(-25,000, -10,000) 


0.17±0.15 


-1.49±0.19 


151.9°H;g 


41 


11 


1.8 10~ 5 


Si II 385.9 


(-25,000, -10,000) 


1.75±4.28 


-5.42±15.1 


140.2°^ 1 8 8 


30 


5 


1.1 10" 5 


Mg II 448.1 


(-25,000, -10,000) 


-0.01±0.12 


-1.13±0.18 


155.8°1|? 


7 


7 


0.34 


I 777.4 


(-25,000, -13,000) 


0.72±0.66 


35.4±293.8 


-89.1 


19 


11 


0.05 


Si III 456.5 


(-15,000, -5,000) 


0.09±0.46 


-1.22±1.02 


154.6°T^ 9 6 6 


11 


4 


0.02 


S II 562.3,544.8 


(-15,000, -9,500) 


-0.24±0.18 


-0.96±0.45 


158.1°!^ 


9 


6 


0.15 


Ca II IR Triplet 


(-24,000, -7,000) a 


0.29±3.6 


-8.5±10.9 


138.3°^ 5 9 ° 


38 


17 


0.0 04 


continuum 


(640, 700) c 


0.11±0.04 


1.1 ± 0.3 


0.0°!:?.! 


46 


23 


0.05 



"Velocity calculated with respect to the restframe wavelength of 854.2 nm. 

h The errors on P. A. are highly non-Gaussian when the slopes of the linear fits are not well constrained. The 
P.A. is periodic with period of 90°. The angle is practically undetermined for Si II 385.9 nm, O I 777.4 nm, 
and the Ca II IR Triplet. 

c Wavelength range in nm 
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Fig. 1. — Spectroscopy of SN 2004dt de-redshifted assuming a recession velocity of 5915 
km s _1 of the host galaxy. The four panels give the total flux spectrum, the polarization, P, 
the reduced Stokes vectors, Q and U, and the polarization angle 6. Typical spectral features 
are marked on the top panel. The lines are produced by material moving at velocities of 
17,000 km s _1 (marked by solid lines) and 11,000 km s" 1 (marked by dotted lines; see Fig. 
3 for a more complete description) in the rest frame of the supernova. Lines from different 
elements are marked with different colors (gray scales). The panels for polarization data 
are binned to 25 A. The wavelength regions used to determine interstellar polarization are 
marked by horizontal lines on the panels for Q and U. Stronger lines tend to be associated 
with larger degrees of polarization except for the O I 777.4 nm line. 
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Fig. 2. — Spectropolarimetry of SN 2004dt displayed in the plane of the reduced Stokes 
vectors Q and U. Note the existence of a dominant axis represented by the solid line. The 
1-a location of the polarization due to interstellar dust is shown as an open circle. The data 
points are binned to 25 A. 
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Fig. 3. — Line identifications in the observed spectrum using the SYNOW code. The only 
ions contributing to the model spectrum (solid line) are Si II, Mg II, O I, and S II. Other 
relevant fit parameters are: The ions are distributed in a shell from 16,000 km s" 1 to 25,000 
km s _1 , with the exception of the S II which is distributed from 11,000 km s -1 to 18,000 
km s _1 ; the location of the photosphere is assumed to be at 11,000 km s _1 with the tempera- 
ture fixed to 11,000 K; the optical depth is assumed to decrease outwards with an exponential 
law with an e-folding velocity of 5 km s _1 ; and the optical depths used in SYNOW are 2.0, 
0.5, 0.5, and 1.0 for Si II, Mg II, S II, and O I, respectively. The unusual appearance of the 
observed Si II profiles are likely due to a combination of the high velocity and overabundance 
of Si II. 
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Fig. 4. — The top panel gives the total flux spectra; the next two panels give the polarization 
spectrum projected onto the dominant axis and the orthogonal axis, respectively. The solid 
vertical lines mark the positions of the lines of Si II, Mg II, Ca II, and O I which are typically 
blue shifted by 17,000 km s . The vertical dash lines mark the positions of the lines of S II, 
Si III, and Fe III, at a blue shift of 11,000 km s" 1 . All the polarized features except O I line 
share the same dominant axis, suggesting that the geometrical structures of Si II, Mg II, and 
Ca II are the same. No polarization is detected across the O I line. Weaker lines of Si III, 
and S II show also polarization along the dominant axis. The data are binned to 25 A. Note 
that the polarization features at 747.0 nm may not be reliable because of telluric absorption. 
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Fig. 5. — The absorption features in the total flux spectrum and the polarization profiles 
along the dominant axis are shown as a function of velocity for the lines of Si II 635.5 nm 
(a, top left), Mg II 447.1 nm (b, top right), O I 777.4 nm (c, lower right), and Si II 385.9 nm 
(d, lower left). The dotted lines show the P-Cygni profiles. The solid histograms shows the 
spectropolarimetry profiles. Note that the polarization is detected in all lines except that 
of O I. Note also that the emission component of the P-Cygni profiles are all blueshifted by 
about 4,000 km s -1 . 
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Fig. 6. — Q-U diagrams of the strongest lines and the corresponding linear fits. The dashed 
lines are the results of the linear fits. The ranges of the data points used for the fits are 
given in Table 2 and are from -25,000 km s -1 to -10,000 km s _1 for these lines. Clockwise 
from top left, the panels are for the lines of Si II 635.5 nm (a, top left), Mg II 448.1 nm (b, 
top right), O I 777.4 nm (c, lower right), and Si II 385.9 nm (c, lower left). Only the Mg II 
line is consistent with a straight line, indicating a simple axially-symmetric geometry with 
no detectable clumping (see text). 
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Fig. 7.— Comparisons of spectra of SN 1997bp, SN 2002bo, SN 2004dt, SN 1990N, 
SN 1994D, and SN 2001el before optical maximum. The data of SN 1997bp were obtained 
at the 2.1 meter telescope of McDonald Observatory through our program of supernova po- 
larimetry (Wang & Wheeler 1997). The data on SN 2002bo and SN 2001el were obtained by 
our supernova polarimetry program at the VLT (Wang et al. 2003a). The data of SN 1994D 
is from Patat et al. (1996). The data of SN 1990N is from Jeffery et al. (1992). The Si II 
635.5 nm lines of SN 1997bp, SN 2002bo, and SN 2004dt are significantly broader than those 
of SN 1990N, SN 1994D, and SN 2001el around optical maximum. The vertical line marks 
the absorption minimum of the Si II 635.5nm line of SN 2001el, SN 1994D, and SN 1990N. 
The Si II 635.5nm lines of SN 1997bp, SN 2002bo, SN 2004dt show significantly larger blue 
shift than those of SN 2001el, SN 1994D, and SN 1990N. 



